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A preliminary investigation of a large class of Hidden Valley models is presented. These models 
are more challenging than those considered in arXiv:0712.2041 [hep-ph]; although they produce a 
new light resonance which decays to heavy standard model fermions, they exhibit no light dilepton 
resonance. A heavy Z' decaying to v-hadrons, which in turn decay mainly to bottom quarks and 
tau leptons, is considered; six case studies are investigated, using a new Monte Carlo simulation 
package. It is found that the one-to-one correspondence of jets and partons is badly broken, and the 
high-multiplicity heavy-flavor signal probably cannot be isolated by counting jets, with or without 
heavy-flavor tags. Instead, other measures, such as counting and correlating vertices or displaced 
tracks, and possibly counting of (non-isolated) muons and use of event-shape variables, should be 
combined with scalar transverse energy and/or missing transverse energy to reduce backgrounds. 
Within the resulting sample, searches for the v-pion mass resonance in both di-jet and single-jet 
invariant mass can help confirm a signal. The best observable in a perfect calorimeter seems to be 
single-jet invariant mass for jets of larger radius (R=0.7), although this needs further study in a 
realistic setting. A more detailed signal-to-background study is needed as a next step, but will face 
the difficulty of estimating the various high-multiplicity backgrounds. 



I. INTRODUCTION 

The "Hidden Valley" scenario [1 , if realized in nature, 
may result in unusual and little-studied phenomena at 
the LHC. In this scenario, the standard model is accom- 
panied by a nearly hidden sector containing light parti- 
cles (a "hidden valley" ) . These particles cannot currently 
be abundantly produced, due typically to an energetic 
barrier or a weak coupling. The increased energy of the 
LHC may greatly enhance their production. The same 
barrier can be traversed in the opposite direction to al- 
low some of these particles to decay visibly to standard 
model particles. A schematic illustration of such models 
is shown in Fig. [I] Examples of hidden valley models 
include the original illustrative classes of models given in 
[T] , along with quirk and squirk models jTJ [^1 El S] and 
a wide class of "unparticle" models [5 a with an added 
mass gap (e.g. [HI E])-. whose signals were discussed in 
detail in [8] . Another related class of examples was stud- 
ied in [TO] E] • Motivation for such sectors is provided 
by, among other possibilities, supcrsymmetry-breaking 
models, which often introduce one or more hidden gauge 
groups. While these gauge groups are normally imag- 
ined to be unimportant at LHC energies, this reflects a 
theoretical bias. Such sectors might also be responsible 
for dark matter, and may have an important role to play 
in other aspects of particle physics [T2], astrophysics and 
cosmology [T5] . 

The main interest of these models for the LHC era is 
that their signatures are often distinctive, and can dif- 
fer from the many supersymmetry, little Higgs, extra- 
dimensional and technicolor signatures that have been 
so often discussed. These include high-multiplicity events 
(generally non-thermal and non-spherical), possibly with 
large missing energy, and exhibiting large event-to-event 
fluctuations. New decay modes for Higgs bosons [TJ |5J 
H4] , supersymmetric particles [H [15] , and top quarks [8] 
often arise. Light neutral resonances are common, lighter 



perhaps than 100 GeV or even 10 GeV. 

A common hidden valley signature, but one which I 
will not address in this paper, is displaced vertices. In 
certain regions of parameter space, some of the new light 
particles have long lifetimes, decaying at macroscopic dis- 
tances. There is no standard model background to com- 
pute or estimate. Key issues associated with such sig- 
nals involve experimental challenges: triggering, detec- 
tor noise, beam halo, pion collisions with detector mate- 
rial, vertex reconstruction, etc. All of these are detector- 
specific, and any study of this signature requires a full 
detector simulation. 

However, it may happen that all the new particles de- 
cay promptly, or that some decay promptly and all others 
are stable and invisible. In this case, identifying the re- 
sulting high-multiplicity and often low-rate signal, over 
a large standard model background, becomes a challenge 
that can be addressed in part through theory and simu- 
lation. Below I will consider models of this type. 

In some classes of hidden valley models [TJ |8] , there is 
a new and frequently-produced particle that often decays 
to electron and muon pairs. It is then relatively easy to 
discover the signal, as emphasized in |16j . Simple and 
rather crude event-shape cuts that remove the largest 
backgrounds may enhance signal-to-background to the 
point that, using the excellent low-px dilepton mass res- 
olution of the LHC experiments, a resonant peak could 
be detected. 

In this paper, I will consider a much more difficult situ- 
ation. I will examine a class of hidden valley models with 
prompt decays and heavy-flavor final states, and with no 
dilepton resonances. Some of these models also have large 
missing energy. All have large event-to-event fluctuations 
in the multiplicity of standard model partons in the final 
state. The background to such signals is difficult to esti- 
mate, because it consists of a cocktail of many different 
processes, none of which can be calculated beyond lead- 
ing order in a s , and few of which can be identified and 
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measured using the data itself. Because of this, it will be 
a considerable challenge to carry out a concrete signal- 
to-background study. The goal of this article is to lay 
the groundwork for such a study, and suggest features of 
the signal which could be used in any search for a hidden 
valley of this type. 

First I will outline the specific model that is chosen as 
an exemplar from this class, and will describe the Monte 
Carlo event generation package used to study it. Then I 
will describe the case studies, examining the basic phe- 
nomenological features of the signal, in Sec. [IT] After ex- 
ploring the basic underlying phenomena, I will consider 
how jets are constructed in the signal. Finding that jets 
are not fully sufficient for interpreting or isolating the 
signal, I will consider other non-standard methods for 
reducing background. Finally, in Sec. |III| I will examine 
the question of how to identify the resonance whose ob- 
servation would confirm the signal, considering both dijct 
and single jet invariant mass. A summary of results and 
some additional comments are given in the conclusion; 
two appendices fill in some details on secondary muons 
and on jet algorithms. 



electroweak doublet and singlet Higgs bosons which mix 
together. These models have been discussed widely (see 
[7] and references therein) but their potential for high- 
multiplicity heavy-flavor final states was only recently 
recognized [T3] |T7] . A third class can include strongly- 
interacting hidden valleys which couple to the standard 
model mainly through the Higgs boson; these have not 
yet been explored fully. 

In this paper I will consider the theory in [T] , as a very 
simple example from the first category. This is a hidden 
valley which closely resembles QCD. To make this study 
especially straightforward, I have chosen a hidden valley 
sector ("v-sector") that, like QCD, has an SU(3) gauge 
group and two light "v-quarks" U and D, with masses ad- 
justed so that the light "v-hadron" mass ratios are those 
of QCD. It is important to emphasize that this model is a 
stand-in for a much larger class of models. Indeed there 
is no reason for the physics of a hidden valley to closely 
resemble QCD, any more than technicolor models should 
closely resemble QCD. However, for initial studies of v- 
sector phenomenology, the case of a QCD-like v-sector is 
simplest to investigate first. This is because the physics 
is easy to understand, and a Monte Carlo event simulator 
is easily constructed. 

As in pQ, where more details are given, I will consider 
such a v-sector coupled to the standard model through a 
broken U(l) gauge symmetry, under which both standard 
model particles and the v-quarks carry a charge. The 
Z' gauge boson of the U(l) will serve to mediate both 
production and decay of particles in the v-sector. 




unstable 
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FIG. 1: A schematic illustration of models in the hidden- 
valley scenario. With sufficient energy, available at the LHC 
but not at LEP, a barrier may be traversed that allows pro- 
duction of new light states in a hidden sector. Dynamics in 
the hidden sector may produce large numbers of particles. 
Some of these new particles may decay back to the standard 
model, often with long lifetimes. 



A. The Models 

Several large classes of hidden valley models share the 
phenomenology of high-multiplicity final states, rich in 
heavy flavor and possibly missing energy. These include a 
wide variety of confining hidden sectors whose light stable 
hadrons are all pseudoscalar and/or scalar mesons with 
comparable masses; an example was given in [TJ. An- 
other class involves weakly-coupled models with multiple 
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FIG. 2: The spectrum of the v-sector considered in this paper. 
In analogy to QCD, all v-hadrons rapidly decay down to v- 
pions and v-nucleons; then the 7r° (and in the B cases, also 
the 7r^ and xjf ) decay more slowly to standard model fermion 
pairs, preferentially to heavy flavor. 

The long-lived v-hadrons of a QCD-like v-sector with 
two light v-quarks U and D are three light v-pions and a 
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heavier v-nuclcon doublet, as shown in Fig. [2} All other 
v-hadrons (such as the v-rho and v-Delta) decay imme- 
diately to v-pions and v-nucleons. For simplicity, it is 
assumed that v-baryon number is conserved, so the v- 
nucleons are stable and invisible. The three v-pions ir v , 
a triplet under v-isospin, consist of a v- flavor off-diagonal 
pair with quantum numbers of UD and DU, analogous 
to the ir^~ of QCD, and a third, the v- flavor diagonal v- 
pion with quantum numbers of UU — DD, analogous to 
the vr°. 

A point of notation: it is natural to name the v-pions 
as (71^, 7r"), in analogy to QCD's pions (Tr*,^ ). This 
notation was used in [T]. However, the use of the ± su- 
perscript proves confusing, because all the v-pions are 
electrically neutral — after all, they are part of a hidden 
sector. To avoid any confusion of this type (and at the 
expense of introducing another), I will call the UU — DD 
state 7T°, but call the UD state ir^, and the conjugate 
DU state 7r,Y. 




FIG. 4: A Z' decays to two v-quarks, which emit v-gluons in 
a v-parton shower. These then are confined into v-pions and 
v-nucleons. Some of the v-hadrons (shown dotted) are stable 
and invisible, but others are metastable and decay, mainly to 
bb. 




FIG. 3: The 7rJJ decays via a Z' to heavy flavor. The 7r^ v , 
if unstable, decays through a v-flavor-changing interaction to 
the same final state. 

If the third component of v-isospin Zjj is conserved, 
then the 7r^ is stable and invisible, but the breaking of 
total v-isospin allows the 7r° to decay via a Z' back to 
standard model particles, as shown in Fig. [3j Helicity 
suppression assures the spin-zero 7r° decays mainly to 
heavy fermions (for the same reason that 7r + — > [i^v 
decays dominate over 7r + — * e + v in QCD); branching 
fractions are roughly proportional to squares of fermion 
masses. In the particular model of [1 , and for light v- 
pion masses, the width of the 7rJ) is 

which has a very strong dependence on model parame- 
ters; here f Vv is the v-pion decay constant, while mz> 
and g' are the Z' mass and coupling. 

It is also possible that the third component of v-isospin 
1% is violated. In this case even the tt£ and can de- 
cay, with widths that are smaller than that of the 7r° by 
a factor which is a dimcnsionless measure of 1% breaking. 
In this article, I simply assume that either (A) 1% is con- 
served (so that the 7r^ is stable and invisible) or (B) 1% 
is badly violated (so that the tt£ decays promptly.) The 
case studies will be divided into "A cases" and "B cases" 
according to this distinction. 

The basic production process for these particles is 
shown in Fig. ffl It involves qq — > Z' — > QQ, where Q is a 



v-quark. The v-quarks undergo a parton shower through 
v-gluon emission, following which they are confined by 
strong v-intcractions into v-hadrons. These v-hadrons 
decay down to v-nucleons and v-pions, and some of the 
v-pions may then decay visibly to standard model parti- 
cles. 



Cross Section in fb 



1000 
100 






10 

1 


— = 20 7>V\>, 

g' 




0.1 







1 2 3 4 5 
Z' mass 



FIG. 5: Cross-section for v-particle production via a Z' in 
the QCD-like model of [I] . The cross-section in other models 
may easily differ by an order of magnitude; see text. 

In the study below, I will consider a Z' with mass 3.2 
TeV. The v-pion masses will range between 50 and 200 
GeV. In this case the v-pions decay promptly, and the 
production cross-section is expected to be of order 10- 
100 fb in the model of pQ, see Fig. [5j In other models the 
cross-section could be different by a factor of 10 or so, 
larger or smaller, due for example to different Z 1 charge 
assignments, or to a different number of colors in the 
v-sector. 

In summary, the model considered below has gauge 
group [SU(3) x SU(2) x U(1)] S m x U(l)' x SU(3) V , with 
the U(l)' broken at the few TeV scale, the SU(3) V group 
confining at the few hundred GeV scale, and two light 
v-flavors of v-quarks U and D. Standard model fermions 
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and v-quarks all carry some charge under the U(l)', al- 
lowing the Z' to serve as a communicator between the 
two sectors. 



— scaling up the masses and momenta of the 
final-state QCD hadrons (with 7r°'s undc- 
cayed) by the factor R, and renaming them 
as v-hadrons. 



B. The Hidden Valley Monte Carlo 0.5 

The event simulator HVMC 0.5, upon which all stud- 
ies in this paper are based, is described in this section. 
(A more general Monte Carlo simulator has been devel- 
oped with S. Mrenna and P. Skands [13], and studies 
based upon it will be presented elsewhere.) HVMC 0.5 
[19] is built on existing tools, which are rather easy to 
modify for current purposes. In particular, elements of 
PYTHIA [20 are strung together to simulate a v-sector 
which is isomorphic to three-color two-flavor QCD, with 
all masses and other dimensional quantities scaled up, 
relative to QCD, by a constant factor R. 



Rescale QQ energy 
? down by factor R 








v-pi-zero — 





FIG. 6: The algorithm used in HVMC 0.5. 



• The decay of the v-hadrons to standard model par- 
tons is simulated using PYTHIA decay routines. 

• The decays, showering and hadronization of the 
standard model partons and the simulation of the 
underlying event proceed using the usual PYTHIA 
routines. 

The resulting final states consist of standard model 
hadrons, photons and leptons, along with stable neutral 
v-hadrons that escape undetected. All results in this pa- 
per are based on analysis of the final state hadrons, pho- 
tons and leptons without accounting for detector effects, 
other than geometric acceptance, except where otherwise 
noted. 

For the sake of clarity (though the actual effect on the 
studies below is small) it should be noted that two-flavor 
QCD as simulated in this way is not quite a consistent 
model. The tuning of PYTHIA to match existing data 
on hadronization, branching fractions, etc., is not correct 
for a two-flavor model. The iso-singlet would-be Nambu- 
Goldstone boson rj v is now affected by the anomaly and 
takes the place of the rf . There are small effects on the 
nucleon mass from the slight differences in the running 
coupling that are similarly ignored. But these issues are 
of minor impact on the phenomenology and of minor con- 
cern for the current studies. As there is no reason to ex- 
pect the hidden sector in nature to be of exactly the form 
considered here, the aim of this paper is not precision but 
rather phenomcnological and experimental guidance, in 
search of robust analysis strategies. 



A v-sector with three colors, two flavors and confine- 
ment scale A„ has v-pions with mass m nv — m^R, where 
R = A v /Aqcd, and Aqcd is the QCD confinement scale. 
It also has nucleons with similarly scaled-up masses. The 
n' v (the iso-singlet pseudoscalar of 2-flavor QCD) has its 
mass set to m v R. Then, given a mass M for the Z', the 
simulation of events proceeds as in Fig. [6] 

• The process qq — > Z' — > QQ is simulated, where q 
is an ordinary quark and Q is a v-quark, using the 
PYTHIA routine for qq -> Z' -> //. 

• The v-parton showering and v-hadronization of the 
QQ system is simulated. This is done by 

— scaling down the energy of the QQ system 
from its original energy Eq to the energy E = 
E /R; 

— simulating QCD parton showering and 
hadronization (with the number of light fla- 
vors set to 2) of an ordinary quark-antiquark 
system with center-of-mass energy equal to 
E; 



II. THE CASE STUDIES 
A. Preliminaries 

The studies below will all involve decays of a Z' of mass 
3.2 TeV to a hidden valley sector. In many models the 
Z' will already been discovered in its decays to dilepton 
final states. However, knowledge of its presence and of 
its mass do not significantly aid in uncovering the hidden 
valley signal, because of the latter 's complexity. In other 
models, the branching fraction of the Z' to dileptons will 
be too small, and the Z' will not yet have been identified 
when the hidden valley signal is sought. 

As noted in Sec. |I A[ signal cross-sections in the 10-100 
fb range are consistent with LEPI and LEPII constraints. 
The cross-section is easily changed by an additional fac- 
tor of 10, without altering the observable phenomenology 
in any other way, by adjusting the U(l)' coupling con- 
stant g' (see Fig. [BJ . For simplicity, I will study samples 
of 1000 events, such as might be obtained in a real LHC 
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FIG. 7: A schematic view of a typical event from case Al. 
The view is along the beampipe. Charged tracks are shown 
in the "tracker", the central disk in the figure. Since there is 
no magnetic field in this event display, all tracks with pr < 3 
GeV have been removed; grey-level corresponds to pr, with 
hardest tracks shown in black. Neutral hadrons and photons 
are indicated as outward-pointing lines starting at the outer 
edge of the tracker, and calorimeter energy in azimuthal an- 
gular bins of width 2n/60 are shown as bars at the outer edge 
of the tracker. 



analysis. All plots shown below, except where noted, show 
results for 1000 simulated signal events. 

The case studies are distinguished by the masses of 
the v-pions and by whether the ir£ are stable or decay 
promptly. In Table |T] the cases are listed. In addition 
to the masses and decay settings, the table shows the 
average multiplicity of visibly-decaying v-pions and some 
kinematic information: the transverse calorimeter energy 
Ht, the invariant mass M 4 of the four highest-p^ jets, 
and the average missing transverse momentum (MET, or 
$t)- The quantities Ht and Ifix are computed here using 
scalar and vectorial sums of the pr of all calorimeter 



towers, 

H T = ]T \p T \e(\p T \-5GeV) 0(H-3) ; (2) 



E 



Pt 



(3) 



The calorimeter towers combine the 3-momenta of the 
various hadrons, electrons, photons and muons in 0.1 x 
0.1 bins in pseudorapidity r\ and azimuthal angle </>. For 
Ht I include only towers with pt > 5 GeV and |7j| < 3, as 
indicated by the functions, to reduce significantly the 
impact of the underlying event. (However, for a fully re- 
alistic study, an Ht variable built from the reconstructed 
jets might be much more robust; in these models, a vari- 
able H^ et \ defined as the scalar summed pt of all jets 
with pt > 25 GeV and \tj\ < 3, takes values about 10 
percent larger than Ht defined above.) Meanwhile M4, 
the invariant mass of the four highest-py central jets, is 
built from jets defined using the midpoint cone algorithm 



of cone radius 0.4; see Sec. IID and Appendix |B| below 




Case 


7T„ 




R 




Ht 


M 4 






stable? 


(GeV) 




decays 


(GeV) 


(GeV) 


(GeV) 


Al 


Yes 


50 


368 


4.0 


667 


590 


318 


A2 


Yes 


120 


883 


2.4 


765 


667 


400 


A3 


Yes 


200 


1470 


1.5 


886 


770 


459 


Bl 


No 


50 


368 


10.3 


1650 


1427 


214 


B2 


No 


120 


883 


6.1 


1835 


1562 


182 


B3 


No 


200 


1470 


3.9 


2248 


1810 


145 



TABLE I: The case studies, showing the stability of the , 
the mass of the v-pion, the ratio R = A v / Aqcd , the average 
number of visible v-pion deays, and the average Ht, Mi, and 



FIG. 8: As in Fig. [7] a schematic view of a typical event from 
case Bl. 

Note the obvious progressions in the table. Comparing 
Al, A2 and A3, one sees the decrease in the multiplicity 
of v-pions; the same trend appears in Bl, B2 and B3. 
Meanwhile the B cases, with a decaying ir£, have roughly 
triple the number of visibly decaying pious, much higher 
visible energy, and much less compared to the A cases 
for the same v-pion mass. For illustration, event displays 
of one event each from cases Al, Bl and A3 are given 
in Figs. [7||9] However the reader should bear in mind 
that event-to-event fluctuations in appearance are much 



G 



larger here than in most standard model backgrounds or 
tradition al n ew signals such as gluino production. 




FIG. 9: As in Fig. [7] a schematic view of a typical event from 
case A3. 
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Missing Energy vs Scalar Transverse Calorimeter Energy 



FIG. 10: The distribution of missing transverse momentum 
(MET) versus Ht in TeV. These quantities are defined in 
Eqs. j2j-(|j. 



Figure 
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shows Ht versus $t for the various cases. 
The difference between the A cases, where roughly 2/3 of 
the v-pions are stable and escape undetected, and the B 
cases, where the v-pions all decay promptly and most of 
the debris from the Z' is observed, is obvious. Occasional 
v-baryons (stable and invisible in all present case studies) 
can provide some fx even in the B cases. There is addi- 
tional and sometimes substantial $t from secondary neu- 
trinos produced in semileptonic decays of b and c quarks 
and especially in t decays. 

Already from these plots, one sees clearly that the A 
cases will have much larger standard model backgrounds 
than the B cases. The B cases are more similar to 
those studied in |16J . though with a higher invariant 
mass, lower rate, and no dilepton resonance. I will 
show in Sec. II D that they have many reconstructed 
hard jets. Backgrounds are high- multiplicity QCD events 
with many b quarks, including ttbb, bbbbjj, ttti, ttZ, etc. 
These signals tend to be in the few pb range or less, 
and will be greatly reduced by an Hx cut at, say, 800 
GeV. The A cases, by contrast, despite their large fx, 
are often in the same overall kinematic regime as rela- 
tively low-energy standard model processes with a few 
jets and fx- These much larger backgrounds include W 
or Z plus jets (especially heavy flavor), tt plus jets, ttW 
or tiZ, tibb, etc. Cuts on H T and fix that have high effi- 
ciency for the signal will still leave considerable amounts 
of background behind. 

In contrast to standard model backgrounds, triggering 
should not be a problem for either the A or B cases. Most 
of the events in this signal will pass the various jet(s) or 
jet(s)-plus-^r triggers, with the latter being most effi- 
cient for the A cases. Even in the A cases, due to sec- 
ondary muons, many events will also pass dimuon and 
muon-plus-^T triggers; see Appendix [A] The trigger will 
mainly remove unspectacular, low-visible-energy events, 
which are rare in the B cases and consist of a large mi- 
nority in the A cases. But the events that fail these trig- 
gers are precisely those which would be especially difficult 
to distinguish from standard model background off-line. 
Conversely, the events which are most distinctive — with 
multiple acoplanar high-p<r jets and possibly large fx ~ 
will be among those which will pass the trigger. For this 
reason, the effect of the trigger is likely to be relatively 
minor, compared to the other issues addressed below. 



B. The v-pions 

In these case studies, the Z' decay produces a sub- 
stantial number of v-pions, organized into two rather fat 
v-jets. All of these v-pions decay visibly in the B cases, 
while about a third are visible in the A cases. (Actually 
the fraction is a bit larger in these QCD-like models, due 
to v-isospin violation that, as in QCD, biases the decay 
of the r] v toward 7r°s.) The number of visible v-pions is 



7 




10 15 20 25 
tt v-Pions 




10 15 20 25 
tt v-Pions 




10 15 20 25 
tt v-Pions 



10 15 20 25 
tt v-Pions 




10 15 20 25 
tt v-Pions 



10 15 20 25 
tt v-Pions 




Number Distribution of Visibly -Decaying v-Pions 



p_T (TeV) p_T iTeV) 

Visibly -Decaying v-Pion Transverse Momentum Distribution 



FIG. 11: For the case studies, with 1000 events, the distribu- 
tion of the number of visibly-decaying v-pions, each decaying 
to two standard model particles. 



FIG. 12: As in the previous figure, the pr distribution of the 
visibly decaying v-pions, in TeV. 



C. v-Pion Decay Products 



shown in Fig. |11| Since the only difference between case 
Bl and case Al is the stability of the 7r^ v , the distribu- 
tion of visible v-pions in case Bl equals the distribution 
of all v-pions, visible and invisible, in case Al. The same 
applies for B2 and A2, and for B3 and A3. 

In Fig.[l2]are shown the pt distributions of the v-pions; 
notice few have pt < 50 GeV, and most are relativistic. 
Fig.[l3]shows the pr distribution of the highest-py visibly 
decaying v-pion, which is almost always relativistic, often 
with a boost factor above 3. This is relevant because 
the majority of events have a highly-boosted v-pion whose 
decay products are separated in r\ and <fi by AR < 0.5. ; 
here AR = \J (A77) 2 + (Ac/)) 2 as usual. We will see in 
Sec. |II D| and Sec. |III| that the decay products of high-p^ 
v-pions often are merged into single jets. 

As an aside, let us note that the multiplicity distri- 
bution and pt distribution of the v-pions, and the dis- 
tribution of final-state quark and lepton flavors, is quite 
different from a thermal distribution, as one would (at 
least naively) expect in a black hole decay |HJ [52] , an- 
other potential source of high- multiplicity events. Also 
the events are not in general spherical (see Figs. [7jj9]) in 
constrast to expectations for black holes. (See however 
Sec. IV D of 0.) 



As noted in Sec. I A the decay rates of the 7r° (and the 



7r^ v if unstable) to quarks and leptons are roughly pro- 
portional to the square of the final-state fermion masses. 
The fermion mass used should be evaluated at the v-pion 
mass scale. For simplicity the relative branching ratios to 
quarks and leptons are taken to be those of a Higgs boson 
of the same mass as the tt v . Since the tt v does not de- 
cay (at tree level) to WW and ZZ, these decay channels 
are first removed before the branching fractions are com- 
puted. In the mass ranges considered here, the unstable 
v-pions decay mainly to bb pairs with a large branching 
fraction (~ 90% for the lighter v-pions) , with the remain- 
der going mainly to t + t~, cc and gluon pairs. The num- 
ber of final-state (short-distance) standard model par- 
ticles is simply double the number of visibly-decaying 
v-pions. Note from Table |T] that the average multiplic- 
ity of final-state partons ranges from 3 in case A3 to 20 
in case Bl. Also, note in Fig. [TT] the wide fluctuations 
and the long tail, which reaches 42 in case Bl and even 
in case A3 extends to 12. (However, a fraction of these 
partons have low transverse momentum, as will be seen 
below.) The event-to-event fluctuations in the multiplic- 
ity of final-state partons are exceptionally high compared 
to most new-physics signals. This is part of what makes 
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FIG. 13: As in the previous figure, the pr distribution of the 
hardest (i.e., highest-pr) visibly decaying v-pion. 



FIG. 14: For the case studies, over 1000 events, the pr dis- 
tribution of the daughters (mainly b and b quarks) of the 
v-pions. 



this signal challenging. 



1. Bottom Quarks and Taus 

Most of the final state partons, especially for lighter v- 
pions where decays to gluons are suppressed, are b quarks 
and antiquarks. The pt spectrum of the central (\rj\ < 2) 
b quarks is shown in Fig. 14 For lighter v-pions, the 



higher multiplicity of b quarks is somewhat compensated 
by their lower px, which makes them less likely to pro- 
duce jets above kinematic cuts and to decay with de- 
tectable vertices. The fraction of bottom quarks that are 
central (\rj\ < 2) and hard (px > 50 GeV) varies from 
about 55% for cases Al and Bl to nearly 85% for cases 
A3 and B3. 

Taus are produced in roughly ten percent of the v- 
pion decays, and are common in these events. One could 
imagine that central hadronic taus could be useful in 
identifying this signal. However, in the present studies, 
few events have more than two tau leptons, which is not 
enough to be unusual (given tt rates). In addition, one 
or both tau's will sometimes fail isolation requirements, 
either because they are too close together (as we will see 
below) or because of the high multiplicity environment 
in which they are produced. For this reason, the num- 



ber of taus identified is likely to be too small for it to 
play a role in extracting the signal. But it should be 
noted that in other hidden valley models, where the r- 
to-6 ratio might be enhanced, the role of taus in signal 
identification might be more important. 



2. Secondary Muons and Electrons 

In these models, unlike those of [16], electrons and 
muons do not provide a direct handle for discovering the 
v-hadrons directly. The branching fraction of the v-pion 
to muons is tiny (unless the v-pion is lighter than 2nib, 
in which case it will be very long lived). However, be- 
cause these events have high multiplicity, and because 
the v-pions decay mainly to b, c and t, which in turn can 
produce light leptons (generally non-isolated), it is very 





Al 


A2 


A3 


Bl 


B2 


B3 


Fraction of events 
with 3 or more muons 
(p T > 3 GeV, \n\ < 2.5) 


0.09 


0.04 


0.02 


0.46 


0.28 


0.11 



TABLE II: Fraction of events with multiple muons. All cases 
studies have 1000 events. 
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common for one or more electron or muon to be produced 
as a secondary. The presence of these light leptons could 
assist with reducing backgrounds. 

Since the leptons in question are not typically isolated, 
I focus on muons, which are easier to identify. The num- 
bers of events with three or more muons that have pt > 3 
GeV and \r)\ < 2.5, with no isolation requirement, are 
shown in Table |Tij Note the large number in the B cases. 
This implies that requiring multiple muons may serve as 
one of several useful criteria for selecting events for anal- 
ysis. Unfortunately the number of muons is too small to 
be of much use in the A cases. 

Additional plots related to lepton distributions appear 
in Appendix [A] 
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Number Distribution of Jets (p_T>50 GeV, |efa|<3) 

FIG. 15: For the case studies, with 1000 events, the number 
of jets (formed using the midpoint cone algorithm with cone- 
size AR = 0.4; see text for more details) with px > 50 GeV 
and \n\ < 3. The number of jets is considerably smaller than 
twice the number of visibly-decaying v-pions, Fig. [TT] 



D. Jets 

Typically one characterizes events on the basis of "ob- 
jects", where the objects include electrons, muons, pho- 
tons, hadronic taus, and jets, which may be tagged or 
untagged. Since the majority of the many jets in the sig- 
nal are from 6's, one might expect about half of them on 
average to be tagged, with a few events containing an ex- 
ceptional number of tags. One might expect these events 
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Number of Partons and Partonic Jets 
vs. Number of Hadronic Jets 



FIG. 16: For the A cases, with 1000 events, the number of 
partons versus the number of hadronic jets (left plot) and 
the number of partonic jets versus the number of hadronic 
jets (right plot). Jets are formed using the midpoint cone 
algorithm with cone-size A_R = 0.4; see text for more details. 
Cuts of pt > 50 GeV and |?7| < 3 are imposed on both partons 
and jets. One sees that that partonic and hadronic jets are 
in correspondence, but partons do not correspond as well to 
jets. 



to be the ones that stand out above standard model back- 
ground. 

This expectation is not entirely wrong, but it is also 
too naive. In particular, the standard jet-parton corre- 
spondence does not work in this signal. As we will see, 
the number of well-reconstructed and taggable jets is typ- 
ically considerably smaller than the number of b quarks. 
Many jets contain two or more b quarks. While the tag- 
ging efficiency may be somewhat higher in such jets, that 
the number of tagged jets obviously cannot be excep- 
tional if the number of jets itself is not exceptional. 



Instead, here and in Sec. HE I will argue that treating 
jets as objects, characterized as either "tagged" or "un- 
tagged", would throw away crucial information needed to 
separate this signal from background. A substantial frac- 
tion of the jets in this signal are not standard jets, and 
it appears that this fact may be critical in suppressing 
backgrounds. 
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vs. Number of Hadronic Jets 



FIG. 17: As in the previous plot, for the B cases. 



In the plots shown in the main part of this article, 
the midpoint-cone jet algorithm is used, with cone size 
0.4; more details on the parameters chosen are given in 
the appendix. Changing parameters, or choosing other 
algorithms, will change the details, but as argued in Ap- 
pendix [S] will not change the main conclusions of this 
section. The same may not be said for finding the v-pion 
resonance, however; see Sec. |III| 

In Fig. 15 is shown the number distribution of jets 
per event with \rj\ < 3 and pt > 50 GeV. The average 



number of jets is rather large but not spectacular in the 
B cases, and not very large in the A cases, which do 
not even have a substantial tail on the high side. This 
means one cannot find this signal by simply demanding 
large numbers of hard jets; in the A cases, a requirement 
of more than six jets removes most of the signal, and 
preserves only half of the signal in the B cases. Recall 
that these signals are in the 10-100 fb range, whereas 
multi-6 backgrounds with 8 or more jets, from ttbb, tiZ, 
etc., are in the 1-10 pb range. A hard cut on the number 
of jets cannot be afforded. 

Note also that in all cases the average number of jets 
is significantly lower than twice the average number of 
visibly decaying v-pions, shown in Fig. |ll| The typical v- 
pion is not producing two jets. It is important to identify 



the reason for this. 

Let us first quantify the degree to which the jets do not 
correspond well to the partons in the event. (In this sec- 
tion, "partons" refers to v-pion daughters, which appear 
at short distance; it does not refer to partons emerg- 
ing through subsequent showering.) That there is a mis- 
match is hardly surprising, given the cluttered nature of 
these high- multiplicity events. In the left-hand plots of 
Figs. 16 and 17 are shown the number of partons ver- 
sus the number of jets; here pr > 50 GeV and \rj\ < 3 
for both partons and jets. Notice these often differ by 
as much as a factor of 1.5 to 2. This breakdown of the 
jet-parton correspondence is natural and indeed has been 
seen before; it will certainly can occur in tttt events, and 
in events with highly boosted massive particles, such as 
W's, Z's, h's and t's. In this signal, however, it can 
become extreme. 

Fortunately, there remains a close connection between 
the clustering of hadrons and the clustering of partons. 
This is partly due to the fact that the final-state quarks 
are all produced in the decays of the color-singlet v-pions, 
which limits the radiation of gluons at large angles. In the 



right-hand plots of Figs. 16 and 17 is shown the relation 
between jets of hadrons and jets of partons. Here, the 
hadrons are clustered according to an algorithm, the par- 
tons (i.e., the short-distance v-pion daughters) are clus- 
tered according to the same algorithm, and the results 
are compared. Clearly the correspondence of hadronic 
jets and partonic jets is much closer than that of hadronic 
jets and partons themselves. In other words, even in this 
signal, a fixed algorithm applied at the hadron level gives 
jets that do correspond well to the jets obtained by ap- 
plying the algorithm at the parton level. It is shown in 
Appendix B that this result is robust for both cone and 



kx algorithms; compare the figures above with Figs. 36 
and |37| In particular, although cone and kx algorithms 
will find different jets in general, both algorithms find the 
same jets in short- distance partons as they do in hadrons. 

Since the partonic jets are a good surrogate for the 
hadronic jets, one can take a short-cut to learn about 
the failure of the jet-parton correspondence. To deter- 
mine precisely the number of b quarks in each jet, one 
should examine the B mesons within the hadronic jets, 
but this is technically tedious and has subtleties. Instead, 
one can examine the number of partons collected within 
partonic jets in various pt ranges, as shown in Fig. |18| 
This provides sufficient information to illustrate the key 
phcnomcnological points. 

Several effects tend to cause multiple partons to be 
combined into a single jet. First, a high-py jet is likely 
to be a single boosted v-hadron, and so contains two 
partons. This effect becomes substantial for boost factors 
above about 3 or 4. This is visible in the right-hand 
plots of Fig. [18] where it can be seen that jets with p? > 
200 GeV have a substantial probability to contain two 
partons, ranging from 1/3 for the higher-mass v-pions 
of cases A3 and B3 to 3/4 for the lower-mass v-pions of 
cases Al and Bl. 
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FIG. 18: For the case studies, with 1000 events, the number 
of partons inside of partonic jets (which correspond well to 
hadronic jets) for jets with \n\ < 3 and pr in three ranges: 
50 < p T < 100 GeV (left plots), 100 < p T < 200 GeV (middle 
plots) and p T > 200 GeV (right plots). 



Also, the fat v-hadronic jets from the Z' decay tend to 
throw multiple v-hadrons into the same region of rj and 
</>, so the probability that partons from different v-pions 
are nearby in i] and <fi is non-negligible. Furthermore, 
the presence of soft partons from the softer v-pions tends 
to increase the probability that harder partons will be 
merged by a jet algorithm. These combined effects can 
be seen in Fig. [18] which shows that a significant number 
of jets contain 3 or more partons, even as many as 6 or 
so, for the cases with lower v-pion mass and consequent 
higher-multiplicity. Of course the effect is more dramatic 
for the B cases. 

Table |III| provides some information about b quarks 
with px less than 30 GeV, which often cannot generate 
a clean jet and are rarely taggable. Recall that the mul- 
tiplicity distributions in these signals have long tails, so 
the average number is less than half the maximum. This 
completely different effect also tends to reduce the num- 
ber of jets relative to the number of v-pion daughters. 
Particularly in the case of lighter v-pions, these low-px b 
quarks contribute additional sources of confusion for jet 
reconstruction, as well as adding tracks and neutrals in 
the few GeV range but without providing a detectable 
vertex. 

Altogether, this means that the number of jets is sig- 
nificantly less than the number of partons. This is a 
bit disappointing, as the high multiplicity of partons is 
a unique and striking feature of the signal. The mere 
counting of jets, even with heavy-flavor tags, is unlikely 
to be enough to separate signal from background, espe- 
cially in the A cases. 



E. Beyond Jets: Vertexing and Tracking 

To identify this signal, it seems likely that tagging of 
individual jets is not enough. By definition, the num- 
ber of heavy-flavor-tagged jets cannot be larger than the 
number of jets. But the number of B mesons can greatly 
exceed the number of tagged jets, as suggested in Figs. [16] 
and [17] In other words, although these events do not 
have an exceptional number of taggable jets, often four 
or less in the A cases, they do have an unusual number 
of B mesons. Thus to distinguish the signal from back- 
ground, it is essential to detect as many vertices from the 
B mesons as possible. 

More precisely, a remarkable feature of this signal is 
the number of vertices and the distinctive correlations be- 





Al 


A2 


A3 


Bl 


B2 


B3 


Fraction of b quarks 
with p T < 30 GeV 


0.27 


0.11 


0.05 


0.25 


0.10 


0.04 


# b quarks per event 
with p T < 30 GeV 


1.86 


0.44 


0.12 


4.49 


1.02 


0.28 



TABLE III: Average distributions over 1000 events of soft 6 
quarks. 
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tween vertices and hadronic jets. Most of the high-p^ jets 
contain more than one B-meson vertex. This can occur 
in some standard model backgrounds (through boosted h 
decays, boosted Z decays, and most commonly through 
g — > bb splitting) but the probability of having two or 
more jets with multiple B mesons is low, and the proba- 
bility of having additional B mesons in the event is also 
low. (Of course, even a single B may produce a sec- 
ond vertex when its daughter D meson decays, but the 
kinematic correlations between the parent and daughter 
vertex are distinctive and different from those of the two 
B mesons from a tt v decay.) It thus appears that moving 
beyond "tagged jets" and "untagged jets" as the basic 
objects of analysis is important for separating signal and 
background. 
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FIG. 20: As for the previous figure, for di-gluon events 
(left-hand plot) and for it (right-hand plot), generated with 
PYTHIA, with radiation and the underlying event. In both 
cases, the samples have Vs > 1 TeV. 
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FIG. 19: For the case studies, over 1000 events, the number 
of displaced tracks versus the total number of tracks. In this 
plot all tracks have \rj\ < 2, pr > 2 GeV, and displacement in 
three dimensions must exceed 300 microns. 

Meanwhile, vertex/jet correlations are not the only 
non-object-based measure that can be useful. Consider 
Fig. |19| in which distributions of the number of tracks 
with |?7| < 2, pt > 2 GeV, versus the number of such 
tracks with three-dimensional impact parameter > 300 
yitm, are shown. This is a measure of both the number of 
B mesons produced and the fraction of tracks that were 
produced in a B meson decay. In particular, notice that 
the slope (the fraction of tracks that are displaced) is 
somewhat larger than in the right-hand plot of Fig. [20] 



which shows tt produced at Vs > 1 TeV. Thus high- 
multiplicity heavy-flavor events will have many tracks of 
which an unusually large fraction will be displaced. 

The clustering of the displaced tracks may also be a 
useful variable, which I have not yet considered. It would 
be interesting to explore a clustering observable acting 
upon them. 

Presumably, the techniques discussed here would be 
useful for many other possible new signals. The need to 
move beyond "tagged" or "untagged" jets is far more gen- 
eral than this particular class of models. It should apply 
in any signal in which a bb pair is produced by a boosted 
particle, such as a Z or h. (For recent relevant work, 
see [13].) Obviously the number of tracks and the frac- 
tion of tracks displaced are blunt instruments, sensitive 
to any process with long-lived particles, whether &'s or 
something exotic and new. The clustering of tracks and 
vertices, however, will be variable from signal to signal. 
For instance, although in the present signal the number 
of vertices is larger than the number of jets, this inequal- 
ity need not hold. In signals with novel heavier long-lived 
particles, which may decay to multiple jets at a displaced 
point inside the beampipe, a number of jets may share 
the same vertex, and then the number of vertices per jet 
may be smaller than one. For example, were case B3 
altered so that the lifetime of the 7r^ v were a few pi- 



coseconds, and were the dominant decay tt^ 



99 



then 



a single decaying n£ v would make two jets, with many 
displaced tracks, emerging from a single vertex. At the 
other extreme, there are models in which a large num- 
ber of light long-lived states are produced, and these can 
have many vertices. Examples would include cases Al 
and Bl with the v-pion mass reduced to 30 GeV and its 
lifetime extended to a few picoseconds. Then the number 
of vertices could be very large due to a large v-pion multi- 
plicity, one vertex per v-pion at the point of its decay, and 
one vertex for each of the daughter B mesons from the 
v-pion decay. This complex of vertexing issues deserves a 
thorough exploration by the fr-tagging community at the 
LHC detectors, including LHCb. 
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FIG. 21: For the case studies, over 1000 events, the (jet- 
level) cluster mass in the two hemi- cylinders divided along 
the transverse thrust axis. Cuts of pr > 25 GeV and \r)\ < 2 
are imposed on the jets. See cautionary remarks in the text. 
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FIG. 22: The (parton-level) cluster mass distribution for the 
Zbbbb background (er ~ 10 fb, 4200 events shown) and Zbbjj 
background (a ~ 2 pb, 10000 events shown). See text for 
more details and cautionary remarks. 



F. A Comment on Event-Shape Variables 

Here I will briefly explore an event-shape variable for 
the signals and for two backgrounds. In [16] it was sug- 
gested that transverse thrust (defined in the two dimen- 
sional plane transverse to the beamline) is a useful vari- 
able for separating signal and background. Another vari- 
able used was "cluster mass" , obtained by dividing the 



cylindrical detector along the plane perpendicular to the 
transverse thrust axis into two hcmicylinders, and com- 
puting the invariant mass of all activity within the hemi- 
cylinder. To reduce the impact of the underlying event 
and of initial state radiation, only calorimeter cells satis- 
fying certain px and rj cuts were used. 

In [TB] , the presence of a dilepton resonance makes even 
a low signal-to-background ratio acceptable. Here, the 
signal is smaller, and a much better signal-to-background 
ratio is needed for the signal to be confirmed (through 
the methods of the Sec. Ill I. This requires that all cuts 
have high efficiency for the signal. There are important 
multi-jet backgrounds (with and without fx for the A 
and B cases respectively) that can be disregarded in the 
case study of [TBJ, but cannot be ignored here. 

A proper background study needs to account for many 
background processes. Below we consider only two, for 
illustration. These are the ~10 fb Zbbbb process and the 
~ 2 pb Zbbjj process (where j is any non-& quark or 
gluon, at least two jets have px > 200 GeV, and the Z 
decays to neutrinos.) [30] The numbers of events shown 
in the plots are 4200 for Zbbb and 10000 for Zbbjj. 

For these two backgrounds, and presumably other 
multi-jet processes that are the dominant backgrounds 
remaining after simple cuts, it appears the transverse 
thrust variable is not helpful. Neither signal nor back- 
grounds resemble back-to-back di-jets, while neither is 
spherical, so no cut removes a large fraction of the back- 
ground without removing most of the signal. But the sit- 
uation with the cluster mass variable is more promising, 
as shown in Fig. [21] for signal and in Fig. [22] for the two 
backgrounds. The cluster masses are here constructed 
from jets, with px > 25 GeV and \r)\ < 2 to reduce 
sensitivity to the underlying event and initial state radi- 
ation. (This is in contrast to [TB], which computed this 
quantity at parton-level; because of jet merging in the 
present signal, such an approach would not be reliable 
here.) The backgrounds (for which jet-parton correspon- 
dence is more likely to hold) are shown at parton-level; 
although this reduces the cluster mass in some events, 
the effect appears small enough to not affect the general 
conclusions below. 

The larger Zbbjj background is clearly the more seri- 
ous problem. It can be reduced if three b tags are re- 
quired of the events, but at a cost of considerable sig- 
nal in the A cases. The cluster mass distribution of the 
signal for the A cases lies underneath the background, 
and appears not to be useful. By contrast, the B cases 
(which, as in [16] . have few invisible final-state particles) 
are much more forgiving, as they are for many variables; 
the cluster invariant mass moves the signal far from the 
backgrounds shown. A loose cut on this variable, com- 
bined with other selection criteria (such as at least three 
6-tagged jets) should help this signal to stand out. 

As an aside, note that the v-pion mass appears visibly 
in the cluster mass distribution. There is a small but non- 
negligible probability that a one or both hemicylinders 
contains only a single v-pion, so that the cluster mass 
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is just the v-pion mass. However the region in which 
this is easiest to see lies underneath the background. It 
seems unlikely that this fact can assist with identifying 
the signal. 

As noted in [IB], the cluster mass variable is not suf- 
ficiently robust for use in a realistic analysis. Studies 
for the present paper have shown marked dependence, at 
the level of 20 percent or more, on the treatment of the 
underlying event and initial state radiation. The plots 
above should therefore be treated with caution. They 
are useful for characterizing differences between signal 
and background, but a more stable version of this vari- 
able should be used in any experimental analysis. 

The conclusions in this section are thus preliminary. 
More robust event-shape variables should be studied, but 
the cluster mass appears to reduce background in the B 
cases. However, the B cases are already distinctive in 
other ways. This variable may be most useful in case 
B3, where the number of jets is not so extreme as in 
case Bl, but the total invariant mass of the jets in each 
hemisphere is still very large. On the other hand, this 
particular variable may not be so useful in the A cases, 
so a different event-shape variable must be sought. 



III. CLINCHING THE CASE: DETECTING THE 
V-PION RESONANCE 

I have discussed a number of features of the signal 
which make its phenomenology atypical. Despite the 
unusual features of the signal outlined above, they are 
not obviously sufficient to allow for easy separation of 
signal from background, if the signal cross-section is in- 
deed 10-100 fb. The standard model backgrounds are 
large and variegated, consisting of tails of distributions 
from a number of different processes. A complete and 
convincing study will be difficult with present tools. In 
any case, it seems unlikely that the backgrounds can be 
understood well enough from data to allow a counting 
experiment, especially in the A cases, where the number 
of jets, tracks, vertices, etc. is not so large. 

Instead, it seems likely that a different strategy is 
needed. Given the unusual features of the signal dis- 
cussed above, one might apply loose cuts on these fea- 
tures that have high efficiency for the signal. (For in- 
stance, one could require substantial |?t and/or Ht, sev- 
eral jets with at least three tagged, indications of many 
displaced tracks and vertices, etc.) The standard model 
background surviving the cuts would not be calculable or 
easily measured, so the signal-to-background ratio would 
not be well-known. But within this enriched sample, one 
could then search for the key kinematic feature of the 
signal - the v-pion resonance - which if observed would 
confirm that new physics is present. 

Simply plotting dijet invariant masses, where the jets 
are selected at random, cannot reveal the v-pion reso- 
nance. The huge combinatoric background, the fact that 
many jets contain multiple fe-quarks, and relatively poor 
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FIG. 23: For the case studies, over 1000 events, the distribu- 
tion of dijet invariant mass for pairs of thin jets (mj < 0.15pr, 
p T > 25 GeV, \rj\ < 3) that are nearby (AT? < 1.2). 



resolution for jet momentum and energy would eliminate 
any signal. 

Since many v-pions are boosted, they often form a sin- 
gle jet, or dissociate into two nearby jets. I will use these 
facts below in studying both dijet invariant mass rrijj and 
single-jet invariant mass rrij below. Clearly both might 
be used; there is no sharp dividing line between them 
in any case, since any jet algorithm separates the two in 
an arbitrary way. (Approaches that avoid this division 
are under study |25j.) It appears that it is important to 
choose one's jets carefully. 

A proper study of these variables would account for 
the finite resolution in jet energy, momentum and mass. 
There are many issues here, some beyond the scope of 
a theoretical investigation. In this study I will include 
decays, showering and hadronization but will work only 
with a perfect calorimeter - perfect except for its granu- 
larity of 0.1 x 0.1 in rj and <j> and its limited pseudorapid- 
ity coverage — and will show that substantial challenges 
arise even before realistic detector issues are accounted 
for. 
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FIG. 24: As in the previous figure, except that the jets have 
Pt > 100 GeV and pairs are within AR = 0.9 of one another. 
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Invariant Mass vs Transverse Momentum of Jets 



FIG. 25: For jets with px > 20 GeV, the distribution of rnj 
versus px, in TeV. 



A. Dijet masses 

Let us begin with dijet invariant mass. A plot of the in- 
variant masses of all pairs of jets above a certain px cut 
would suffer from an overwhelming combinatoric back- 
ground, because of the high jet multiplicity. Instead, it 
is best to use the fact that high-energy Z' decays provide 
a substantial boost to many of the v-hadrons, as we saw 
in Sec. [il] Most high-p^ jets are either single v-hadrons, 
whose decay products have merged, or they represent a 
single quark produced by the decay of a v-hadron whose 
other decay product will lie close by in AR. The dif- 
ference between these two cases is (on average) that the 
former class of jets will have a single-jet invariant mass 
rrij larger than the latter class. In particular, we will dis- 
tinguish between "thin" and "thick" jets, thin jets being 
those with rrij < 0.15pt, and thick jets being those with 
irij > Q.lbpx- (This terminology has been introduced 
in [25].) By selecting thin jets with high px, and plot- 
ting the dijet invariant mass of pairs of such jets with 
AR. not large, one might hope to significantly reduce the 
combinatoric background. (This same method works for 
finding W bosons in high-energy tt events.) 

In Fig. [23] is shown a plot of dijet masses for cone 
jets (defined as described in Sec. |IID and in Appendix 
|b[) I demand that both jets are thin and have px > 25 



GeV and < 3, and require that the two jets have 



AR < 1.2. In Fig. 24 is shown a similar plot with lower 
statistics but lower background, using p T > 100 GeV and 
AR < 0.9. Either approach is a challenge for cases Al 
and Bl, which is not surprising, since reconstruction of 
a 50 GeV resonance using jets is no easy task. As can be 
seen from the plots, the number of v-pions reconstructed 
is disappointingly low; recall there are thousands in the 
data (see Table [I]) Smearing and mismeasurements, not 
included here, will only make matters worse. It would 
appear that dijet invariant mass is not a particularly good 
variable for reconstructing the v-pion resonance. 



B. Single jet masses 



Now let us turn to single jet invariant mass m». In 



Fig. 25 rrij versus px is shown for all jets in these events 
withT^! < 3. The same is shown in Fig. 26 for the 
highest-pT jet in each event. In both classes of plots, 
a band of jets with invariant mass near to the v-pion 
mass is clearly seen. However, this information is not 
entirely observable. High px QCD jets will often develop 
an invariant mass of order 15 percent of their px just 
from the emission of a moderate- kx gluon that lies out- 
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Invariant Mass vs Transverse Momentum of Hardest Jets 



FIG. 26: For the highest-pr jet in each event, the distribution 
of rrij versus pr, in TeV. 



side the parton-shower approximation used in PYTHIA's 
simulation of jets. Also, jets which are too narrow will 
not have a well-measured invariant mass, for various 
detector-related reasons. One might therefore wisely ex- 
clude from analysis any jet whose mass is less than, say, 
15 or 20 percent of its pr, thus excluding the hardest jets. 
Refinements of this measurement deserve more attention 
than can be given here. However, some preliminary in- 
dications are presented in Fig. |27| where the single jet 
invariant mass of all jets with px > 100 GeV is shown on 
the left, and the same with the additional condition that 
the jet be "thick", rrij > 0.15pt, is shown in the second- 
to-left column. Clearly the additional condition has the 
advantage of removing many ordinary isolated jets, re- 
ducing the QCD continuum that peaks at low mass, and 
allowing the signal to stand out more clearly. This more 
refined measure of invariant mass is shown again in the 
right-hand plots of Fig. [27] for the highest-pr and second- 
highest-pT jet in each event. Interestingly, in case Bl — 
the case with highest multiplicity, where the highest-p^ 
jet often is a merging of more than two partons (Fig. 18 ) 
— the second-highest-pT jet shows the v-pion mass more 
clearly. The correlation between the masses of the two 



highest jets can also be a useful variable; a scatter plot of 
the masses of the two highest-p^ jets (if both are central 
and thick) is shown in Fig. 28 Note that in case Bl the 



high multiplicity tends to spread out the peak, whereas 
in case A3 the peak is almost invisible due to low statis- 
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FIG. 27: For the case studies, over 1000 events, the plots, from 
left to right, show (a) the distribution of single jet invariant 
mass rrij for jets with pr > 100 GeV and \rj\ < 3; (b) the 
same, with the additional condition that the jets be "thick" , 
namely rrij /pr > 0.15; (c) the invariant mass of the highest- 
p T jet (if p T > 100 GeV, \t]\ < 3, and m 3 /p T > 0.15), and (d) 
the same as the previous plot but for the second-highest-pT 
jet (with the same cuts.) 



tics, but for the other cases this plot can help reveal the 
resonance. 

The fact that the single-jet invariant mass is a good ob- 
servable in all 6 cases is a simple consequence of the fact 
that relatively light particles are being produced with a 
large boost, through the decay of a heavy Z' . (Clearly 
the same strategy will not work for v-pions produced in 
decays of lighter particles, such as Higgs bosons [U [14] 
or supersymmetric particles |15j.) There are important 
and little-studied backgrounds from all-hadronic decays 
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FIG. 28: Jet-invariant mass (in GeV) for the highest and 
second-highest pr jets in each event; jets are required to have 
p T > 100 GeV, mj > 0.15p T , and |r/| < 3. 



of boosted W's, Z's, and t's, which are produced with 
an enormous rate. These will swamp any new resonance 
unless events are first selected with the unusual features 
of this signal. Should the v-pion mass lie close to 80-90 
GeV or to 170 GeV, the difficulties will be very much 
greater. We must hope nature does not choose this sce- 
nario, or at least provides a large cross-section in return. 



C. Some Improvement Using Fatter Jets 

We have seen that single jet invariant mass recon- 
structs more v-pions than does dijet invariant mass. By 
widening the jet cone, one might hope to reconstruct even 
more. Here I will explore increasing the cone radius from 
0.4 to 0.7, and will show some increase in efficiency. 



Comparing Fig. 29 to Fig. [25] we see that a larger frac- 
tion of the jets of radius 0.7 are single v-pions, compared 
to those of radius 0.4. This is very clear in the left-hand 
whereas the thickness criterion (the cut 
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plots of Fig. 

in rrij /pt) is essential to remove random jets in the signal 
in Fig. [27] it is less essential, but still effective, with the 
larger cone size. Comparison of Fig. [30]with Fig. [27] and 




Invariant Mass vs Transverse Momentum of Jets 



FIG. 29: As in Fig. [25] but with cone size 0.7 instead 0.4; for 
jets with pr > 20 GeV, the distribution of rrij versus pr, in 
TeV. 



of Fig. 31 with Fig. [28] shows that the larger cone size 
generally allows a marked improvement in both efficiency 
and resolution, for both the hardest and second-hardest 
jet in the event. 

The one interesting exception, among the case studies, 
is case Bl. Here the number of partons in the final state 
is so large that confusion background dominates. Most 
high-pT jets contain multiple b quarks, and the hardest 
jets tend to contain more than 2, a tendency already vis- 
ible in Fig. 
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With a cone size of 0.4, rather few of 
the high-pT jets are single v-pions, and a larger cone-size 
makes this problem worse. It appears it is best in this 
case to work with the jets of radius 0.7 that do not have 
the highest pr, and even then the background from ran- 

fourth 
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dom jets in the signal is rather large (see Fig. 
line, second plot from left). Background from standard 
model processes would be very problematic except for the 
fact that case Bl is also the easiest case to separate from 
the standard model using other methods. This case has 
the highest multiplicity of jets (~ 7) and vertices (~ 20), 
the highest multiplicity of secondary muons (~ 3), the 
most tracks (~ 100) and displaced tracks (~ 50), very 
high Ht (~ 1-6 TeV), and a striking event shape (sum 
of the two cluster masses ~ 1.2 TeV). Perhaps this signal 
can even be identified in a counting experiment, where 
the standard model backgrounds can be estimated from 
the data by looking at event samples that share some but 
not all of these striking features. It is conceivable that 
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FIG. 30: As in Fig. [27] but with cone size 0.7 instead 0.4; 
from left to right, (a) the mass rrij of jets with pr > 100 GeV 
and \rj\ < 3; (b) the same but with rrij /pr > 0.15; (c) rrij for 
the highest-pr jet with jets with pr > 100 GeV, < 3 and 
rrij/pT > 0.15; and (d), the same but for the second-highest- 
Pt jet. 



the v-pion resonance can be better identified with a more 
sophisticated variable than single jet mass, looking more 
carefully at the substructure of the jets. (It is even pos- 
sible that, with so many v-pions per event, and with a 
bit more statistics than available here, the v-pion can be 
discovered through its rare tree-level decay to muon pairs 
or its loop-induced decay to photon pairs.) More gener- 
ally, it is important to study further how best to look for 
resonances in very-high-multiplicity signals, such as case 
Bl. 
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FIG. 31: As in Fig. [28] but with cone size 0.7 instead 0.4; jet- 
invariant mass (in GeV) for the highest and second-highest pr 
jets in each event, where jets are required to have pr > 100 
GeV, rrij > 0.15pr, and |r;| < 3. 



D. Comments 

I have given evidence that single jet mass, using a 
larger cone size than typically used for jets at the LHC, 
is the variable to use in searching for the v-pion reso- 
nance. Why are the fatter jets a better choice? If a sin- 
gle boosted v-pion has a boost factor greater than 6, it 
will typically form a single thin jet. Thus, thick jets arise 
from v-pions with a boost below 6, whose daughters typ- 
ically have an opening angle of order 0.3 or larger. This 
is why, if a thickness criterion is applied, jets of radius 
0.7 have a much higher efficiency for containing a single 
v-pion than jets of radius 0.4. 

These conclusions are somewhat suspect, and the plots 
above unrealistic, because no energy smearing or mag- 
netic field were included in the simulation of the calorime- 
ter. A more serious study is needed, using a more com- 
plete detector simulation. Several remarks are in order. 

It will be of considerable interest to learn the jet-mass 
resolution of the LHC detectors. This information should 
be available in the early data from studies of W bosons 
and t quarks in boosted ti events. It will also be inter- 
esting to see plots of the QCD continuum background to 
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jet-mass measurements. 

Since angular resolution is essential for jet mass mea- 
surements, the curvature of tracks due to the magnetic 
field should be removed. Clearly the use of something like 
"particle flow" — combining angular information from 
the tracker as well as the electromagnetic and hadronic 
calorimeters — should significantly aid in improving the 
resolution on jet mass. 

An additional handle for reducing backgrounds to 
boosted v-pions may lie in detecting the substructure 
within the jet, such as recently considered in [23 . This 
is easiest to do with clustering algorithms such as the 
kx algorithm, and using tracking as well as calorimetry. 
In fact, for the most energetic v-pions, whose daughters 
have the lowest angular separation, the granularity of 
the calorimeter, especially its hadronic component, will 
badly degrade jet-mass resolution. In this case tracking is 
crucial (see for example [26.) By studying the locations 
of the highest-pT tracks, one may be able to significantly 
improve angular resolution on the substructure of the jet, 
and improve the invariant mass measurement. These is- 
sues deserve much more attention and will be explored 
elsewhere [2"5] , 

IV. CONCLUSION 

In this paper I have investigated the phenomenology 
of a large class of hidden valley models, illustrated by six 
case studies within a particular model, but intended to 
apply in a much broader context. The models present 
a single new narrow resonance (more precisely, in the 
particular case studied, a triplet of nearly degenerate 
resonances 7r" , 7r^ , ir^ ) that decays promptly, and pre- 
dominantly to heavy flavor. Importantly, there is no 
measurable rate for decay to or e + e~. In this 

sense this class is largely orthogonal to the large class 
studied in [16 . In production mechanisms (Z' decay) 
that lead to high-multiplicity heavy-flavor final states, 
the phenomenology is very different from that of most 
supersymmetry, technicolor, Randall- Sundrum, or little 
Higgs models that have so often been studied. 

The six case studies were divided into the A cases 
(where the 7r^ v are invisible and stable) and the B cases 
(where all three types of n v decay promptly). The main 
results were the following: 

• Missing energy and total transverse energy are use- 
ful global variables. 

— In the A cases, fx and Ht are of the same 
order, many hundreds of GeV for 50 GeV v- 
pions and closer to a TeV for 200 GeV v-pions. 

— In the B cases, the $t is much smaller and the 
Ht is in the 1.5-2.5 TeV range. 

• Multiplicities of v-hadrons, and their daughter 
standard model partons (mostly bottom quarks), 
are very large. 



— In the A cases, the number of v-hadron daugh- 
ters ranges from an average of 3 with a tail to 
12 (for heavy v-pions) to an average of 8 with 
a tail to 24 (for lighter v-pions). 

— In the B cases, the number of v-hadron daugh- 
ters ranges from an average of 8 with a tail to 
20 (for heavy v-pions) to an average of 20 with 
a tail to 42 (for lighter v-pions). 

• However the multiplicity of jets is somewhat 
smaller, because jets often contain multiple v- 
hadron daughters; this is due to the high boost and 
high concentration of the v-pions. 

— In the A cases, this effect is especially impor- 
tant, because by reducing the number of hard 
jets it leaves the signal with larger standard 
model backgrounds. 

— In the B cases, the effect is more pronounced, 
but the resulting multiplicity of jets is still 
large compared to most standard model pro- 
cesses. 

• In the A cases, counting of heavy-flavor-tagged jets 
appears insufficient for separating signal from back- 
ground, because the multiplicity of jets is too low. 
The situation in the B cases is somewhat more 
promising. 

• The signal is particularly unusual in the numbers 
of vertices, of tracks, and of displaced tracks, and 
the clustering of and correlations among jets, tracks 
and vertices. These features may be useful in sep- 
arating signal from pbackground. 

• Other observables, including event-shape variables 
(particularly one similar to the M c i uster variable of 
|16j ) and numbers of secondary muons, might serve 
as additional tools for event selection, though their 
utility needs more study. They may not be helpful 
in the A cases, and they may not be needed in the 
B cases. 

• If the unusual features of the signal can be used 
to obtain a sample with a reasonable signal-to- 
background ratio, an attempt can be made to find 
the v-pion resonance. 

— Because of combinatorics, a naive approach to 
dijet masses cannot work. 

— Since many v-pions are boosted, it is better 
to consider dijet masses of jets which are close 
in r\ and <f>, or single-jet masses of individual 
jets which are "thick" (have a large mass-to- 
Pt ratio.) 

— Single-jet mass for a relatively fat jet defini- 
tion seems to be an observable with a high 
efficiency for the signal. Here R = 0.7 cone- 
jets were shown to be better than those with 
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R = 0.4, but a thorough study of the best 
jet radius was not carried out, nor were other 
algorithms carefully considered. More work 
remains to optimize this measurement [55] . 

Correlations between the single-jet mass of the 
highest- and second-highest-pT jets may be 
useful in reducing background. 

With extreme multiplicity, as in case Bl, jet 
mass for high-pr jets is less useful, because too 
many of these jets contain more than a sin- 
gle v-pion. An alternative observable is not 
yet known, but also such signals are partic- 
ularly spectacular, with low standard model 
backgrounds, and a discovery claim might not 
require reconstruction of the resonance. 

It is possible but rather unlikely that multi- 
plicities and event rates will be high enough 
to allow reconstruction of the n v resonance 



through the rare decays ir v 
77- 



or ir v 



In all cases, the phenomenological issues that arise are 
somewhat unusual, and need to be explored further in 
more fully realistic studies. Generating relevant and re- 
alistic background samples for these unusual signals rep- 
resents a substantial challenge. 

I should again emphasize that the issues and observ- 
ables discussed in this paper are not limited to these spe- 
cific models, but will apply more broadly. In a number 
of classes of hidden valley models, new resonances that 
decay to heavy flavor are produced in high-multiplicity 
environments, sometimes with a substantial boost. This 
can also happen in models beyond the hidden valley sce- 
nario. 

Conversely, it is important to stress that this particular 
class of models is special, and relatively easy, in that there 
is only one v-hadron resonance to be found. In many 
other models, there will be multiple v-hadrons of differ- 
ent masses, and so the resonance signal just described 
will be spread out among several resonances, making it 
harder to extract. One may still hope, however, for an 
enhancement in jet and dijet masses that will be inconsis- 
tent with Standard Model background. And some mod- 
els are even easier to find, as in the example [T] studied 
phenomcnologically in 16J , because of the presence of 
electron-pair and muon-pair resonances. These may eas- 
ily be seen above background even with a rather impure 
sample of hidden valley events. 

Several other classes of hidden valley models with dis- 
tinct phenomenology from those studied here and in |16j 
remain to be investigated. To examine these in detail re- 
quires a significant extension of the current Monte Carlo 
simulation package. An initial extension is now complete 
[18J. Preliminary results on some models, with their own 
unusual signals, will be presented soon |27j . 
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FIG. 32: The distribution of the number of secondary muons. 



APPENDIX A: 



OBSERVATIONS CONCERNING 
LEPTONS 



Plots of the number of central prompt or semi-prompt 
(appearing within the beampipe) muons, with px > 3 
GeV and \<q\ < 2, are shown in Fig. 32 and their px 



distributions are shown in Fig. 33 
are imposed. 



No isolation criteria 



APPENDIX B: CHECKING THE 
CORRESPONDENCE OF HADRONIC JETS AND 
PARTONIC JETS 

Here I will show more carefully that the partonic and 
hadronic jets, constructed with a reasonable algorithm, 
do correspond, as claimed in the main text. 

The jet algorithms used in this study are shown in Ta- 



ble IV All studies use the mult i- algorithm software Spar- 
tyJet of [28], which includes the Fast Jet kx algorithm 
|29j . In the plots within the main text, the midpoint 
merging version of the cone algorithm was used, with 
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FIG. 33: The pr distribution of the secondary muons. 
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FIG. 34: Correlation between the pr of the highest-pr par- 
tonic jet and the pr of the highest-pr hadronic jet, con- 
structed using the midpoint cone algorithm with R — 0.4 
(see Table [TV]) 



cone size R — 0.4. In this appendix, a cone of R = 0.7 
is also used, along with the kr algorithm with the R pa- 
rameter set to give kj- jets of radius approximately 0.4 
and 0.7. 

First, let us check the correspondance between 
hadronic jets and parton jets. One may confirm that 
the pt of the highcst-pT partonic jet matches with the 
Pt of the highest-p^ hadronic jet; this is shown in Fig. |34| 
The psedorapidities match as well. Nor is this an acci- 
dent limited to the highest-p^ jet; the same figure for the 
second-highest pt jet shows that the strong correlation 
between the hadronic and partonic jets persists. 

Now I turn to the kr algorithm. As with the cone 
algorithm we see, in Figs. [36] and [37] a failure of the cor- 
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FastJet kr Algorithm 


R 
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dcut 
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TABLE IV: The parameters used in the cone and fer algo- 
rithms used in this study. 
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FIG. 35: As in the previous plot, for the second-highest-pT 
partonic and hadronic jets. 



rcspondcnce of partons to hadronic jets, and better agree- 
ment between partonic jets and hadronic jets. Figures [38] 



and 39 further confirm the quantitative correspondence 
of the partonic and hadronic jets. Thus the cone and fc<r 
algorithms both allow reconstruction of the partonic jets 
using hadronic jets. 
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These studies have been repeated for jets with larger 
settings of the cone size or R parameter (0.7 and 1.0 for 
cone jets) . The correspondence of partonic and hadronic 
jets continues to hold firm. 



Number of Partons and Partonic Jets 
vs. Number of Hadronic Jets 



FIG. 36: As in Fig. |16| for the kr algorithm with R parameter 
0.52 (see Table [Tv]) 
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FIG. 37: As in Fig. |17| for the fcy algorithm with R parameter 
0.52 (see Table |Tv|) 
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FIG. 38: As in Fig. |34| for the kr algorithm with R parameter 
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